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Abstract. We study how parallel chip-firing on the complete graph 
K n changes behavior as we vary the total number of chips. Surprisingly, 
the activity of the system, defined as the average number of firings per 
time step, does not increase smoothly in the number of chips; instead 
it remains constant over long intervals, punctuated by sudden jumps. 
In the large n limit we find a "devil's staircase" dependence of activity 
on the number of chips. The proof proceeds by reducing the chip-firing 
dynamics to iteration of a self- map of the circle S 1 , in such a way that the 
activity of the chip-firing state equals the Poincare rotation number of 
the circle map. The stairs of the devil's staircase correspond to periodic 
chip-firing states of small period. 



In this paper we explore a connection between the Poincare rotation num- 
ber of a circle map S 1 — > S 1 and the behavior of a discrete dynamical system 
known variously as parallel chip-firing [4, 5] or the (deterministic) fixed en- 
ergy sandpile [2, 15]. We use this connection to shed light on two intriguing 
features of parallel chip-firing, mode locking and short period attractors. 
Ever since Bagnoli, Cecconi, Flammini, and Vespignani [1] found evidence 
of mode locking and short period attractors in numerical experiments in 
2003, these two phenomena have called out for a mathematical explanation. 

In parallel chip-firing on the complete graph K n , each vertex v € [n] = 
{1, . . . , n} starts with a pile of o~(v) > chips. A vertex with n or more 
chips is unstable, and can fire by sending one chip to each vertex of K n 
(including one chip to itself). The parallel update rule fires all unstable 
vertices simultaneously, yielding a new chip configuration U a given by 
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Here 

r(a) = #{v\a(v) > n} 

is the number of unstable vertices. Write U°a = a, and U t a = U(U t ~ 1 a) 
for t > 1. 

Note that the total number of chips in the system is conserved. In par- 
ticular, only finitely many different states are reachable from the initial 
configuration a, so the sequence (J7*<r)t>o is eventually periodic: there exist 
integers m > 1 and to > such that 

U t+m a = U l a Vt > t - (2) 

The activity of a is the limit 

a(a) = lim — (3) 
<-5>oo nt 

where 

t-i 

a t = J>(i7V) 

s=0 

is the total number of firings performed in the first t updates. By (2), 
the limit in (3) exists and equals ^j-(at+m ~ a t) for any t > to- Since 
< a t < nt, we have < a(a) < 1. 

Following [ I L] , we ask: how does the activity change when chips are added 
to the system? If a n is a chip configuration on K n , write a n + k for the con- 
figuration obtained from a n by adding k chips at each vertex. The function 

s n (k) = a(a n + k) 

is called the activity phase diagram of a n . Surprisingly, for many choices of 
a n , the function s n looks like a staircase, with long intervals of constancy 
punctuated by sudden jumps (Figure 1). This phenomenon is known as 
mode locking: if the system is in a preferred mode, corresponding to a wide 
stair in the staircase, then even a relatively large perturbation in the form 
of adding extra chips will not change the activity. For a general discussion 
of mode locking, see [12]. 

To quantify the idea of mode locking in our setting, suppose we are given 
an infinite family of chip configurations 0-2,03,... with a n defined on K n . 
Suppose that < a n (v) < n — 1 for all v G [n], and that for all < x < 1 

— #{v S [n] I a n (v) < nx} — > Fix) as n — > 00 (4) 
n 

for a continuous function F : [0, 1] —> [0, 1]. Then according to Theorem 8, 
the activity phase diagrams s n , suitably rescaled, converge pointwise to a 
continuous, nondecreasing function s : [0,1] — > [0,1]. Moreover, under a 
mild additional hypothesis, Proposition 10 says that this limiting function s 
is a devil's staircase: it is locally constant on an open dense subset of [0, 1]. 
For each rational number p/q G [0, 1] there is a stair of height p/q, that is, 
an interval of positive length on which s is constant and equal to p/q. 
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Figure 1. The activity phase diagrams s n (k) = a(a n + k), 
for n = 10 (top left), 100 (top right), 1000 (bottom left), and 
10000, where o~ n is given by (5). On the horizontal axis, k 
runs from to n. On the vertical axis, s n runs from to 1. 



Related to mode locking, a second feature observed in simulations of par- 
allel chip-firing is non-ergodicity: in trials performed with random initial 
configurations on the n x n torus, the activity observed in individual tri- 
als differs markedly from the average activity observed over many trials 
[15]. The experiments of [ ] suggested a reason: the chip-firing states in 
locked modes, corresponding to stairs of the devil's staircase, tend to be 
periodic with very small period. We study these short period attractors in 
Theorem 18. Under the same hypotheses that yield a devil's staircase in 
Propositon 10, for each q £ N, at least a constant fraction c q n of the states 
{a n + k}fc =0 have eventual period q. 
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To illustrate these results, consider the chip configuration a n on K n de- 
fined by 



cr n (v) 



n 




v - 1 




+ 




.4- 


2 



1, 



, n. 



(5) 



F(x) 



l 

2' 



I < T < 3 
4 — — 4 

X > " 



(6) 



Here [^J denotes the greatest integer < x. This family of chip configurations 
satisfies (4) with 

'o, 

2:r 

1, ^_ 4 . 

The activity phase diagrams of a n for n = 10, 100, 1000, 10000 are graphed 
in Figure 1. For example, when n = 10 we have 

(a(a 10 + k))l° =0 = (0, 0, 0, 0, 1/3, 1/2, 1/2, 2/3, 1, 1, 1) 

and when n = 100, we have 

(a(crioo + k))l% =(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 1/6, 1/5, 1/5, 1/4, 1/4, 1/4, 2/7, 1/3, 1/3, 1/3, 1/3, 1/3, 1/3, 1/3, 
2/5, 2/5, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, 
1/2, 1/2, 1/2, 3/5, 3/5, 2/3, 2/3, 2/3, 2/3, 2/3, 2/3, 2/3, 5/7, 3/4, 3/4, 3/4, 
4/5, 4/5, 5/6, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1). 



As n grows, the denominators of these rational numbers grow remarkably 
slowly: the largest denominator is 11 for n = 1000, and 13 for n = 10000. 
Moreover, for any fixed n the very smallest denominators occur with greatest 
frequency. For example, when n = 10000, there are 1667 values of k for 
which a(a n + k) = ^, and 714 values of k for which a(a n + k) = ^; but for 
each p = 1, ... ,12 there is just one value of k for which a(a n + k) = j^. In 
Lemma 17, we relate these denominators to the periodicity: if a{a) = p/q 
in lowest terms, then a has eventual period q. 

The remainder of the paper is organized as follows. In section 2 we show 
how to construct, given a chip configuration a on K n , a circle map / : 
S 1 — > S 1 whose rotation number equals the activity of a. This construction 
resembles the one-dimensional particle/barrier model of [10]. In section 3 we 
use the circle map to prove our main results on mode locking, Theorem 8 and 
Proposition 10. Short period attractors are studied in section 4, where we 
show that all states on K n have eventual period at most n (Proposition 16). 
Finally, in Theorem 23, we find a small "window" in which all states have 
eventual period two. 

In parallel chip-firing on a general graph, a vertex is unstable if it has 
at least as many chips as its degree, and fires by sending one chip to each 
neighbor; at each time step, all unstable vertices fire simultaneously. Many 
questions remain concerning parallel chip-firing on graphs other than K n . 
If the underlying graph is a tree [ I] or a cycle [7], then instead of a devil's 
staircase of infinitely many preferred modes, there are just three: activity 
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0, 2 and 1. On the other hand, the numerical experiments of [1] for parallel 
chip-firing on the n x n torus suggest a devil's staircase in the large n limit. 
Our arguments rely quite strongly on the structure of the complete graph, 
whereas the mode locking phenomenon seems to be widespread. It would be 
very interesting to relate parallel chip-firing on other graphs to iteration of 
a circle map (or perhaps a map on a higher-dimensional manifold) in order 
to explain the ubiquity of mode locking. 

2. Construction of the Circle Map 

We first introduce a framework of generalized chip configurations, which 
will encompass chip configurations on K n for all n. To each generalized chip 
configuration we associate a probability measure on the interval [0,2), and 
to each such measure we associate a circle map S 1 — > S 1 . We will define an 
update rule U describing the dynamics on each of these objects, so that the 
following diagram commutes. 

chip configurations on K n — — — > chip configurations on K n 



generalized chip configurations 



measures on [0, 2) 



lifts of circle maps S — > S > lifts of circle maps S S 



generalized chip configurations 



(7) 



measures on [0,2) 



Write C for Lebesgue measure on [0, 1]. A generalized chip configuration 
is a measurable function 

n: [0,1] -> [0,oo). 

Let 

r{rj) = £{x\rj(x) > 1} 

and define the update rule 

I n{x) — 1 + r(rj), rj(x) > 1. 

If a is a chip configuration on K n , we define its associated generalized chip 
configuration ^{a) by 

= a(\nx])+\nx]-nx 
n 

where \y\ denotes the least integer > y. Our first lemma checks that the 
top square of (7) commutes. 
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Lemma 1. U oip = ip oU. 

Proof. Let a be a chip configuration on K n , and write r\ = ip(o~). For v G [n], 
if a(v) > re, then n(x) > 1 for all x G (^-) r]- Conversely, if cr(u) < n — 1, 
then 77(2;) < 1 for all x G (— -, -]. Hence r(r/) = r(a)/n. 
Fix x G [0, 1] and let v = \nx~\ . Then 

£700))(x) = rj(x) - l{,(a)>i} + r(rj) 

a{y) + v — nx r(a) 

= l{cr(t))>n} H 

n n 
J7 + u- nx 

n 

= ip(Ua)(x). □ 



We remark that Lemma 1 would hold also if we defined ip using the 
piecewise constant interpolation . We have chosen the piecewise linear 

interpolation (9) because it is better suited to our construction of the circle 
map, below: the circle map associated to the piecewise linear interpolation 
will always be continuous. 

Next we observe a simple consequence of the update rule (8), which will 
allow us to focus on a subset of generalized chip configurations which we 
call "confined." 

Lemma 2. Let r] be a generalized chip configuration, and let x G [0,1]. 
If r](x) < 2, then 

r(rj) < Un(x) < 1 + r(n). 

Proof. The first inequality is immediate from the definition of Un. For the 
second inequality, if n(x) < 1, then 

Urj(x) = n(x) + r{rf) < 1 + r(n) 

while if 1 < r/(x) < 2, then 

Un(x) = n(x) — 1 + r(i]) < 1 + r(r]). □ 

Definition. A generalized chip configuration rj is preconfined if it satisfies 

(i) r](x) < 2 for all x G [0,1]. 

If, in addition, there exists r G [0, 1] such that 

(ii) r < rj{x) < 1 + r for all x G [0, 1] 
then 77 is confined. 

By Lemma 2, if n is preconfined, then Un is confined. 
Note that from (1) 

Urj(x) = r](x) + r(rj) (modi). 
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Iterating yields the congruence 

U l r]{x) = r/O) + fa (mod 1) (10) 

where 

t-i 

A = 5>(17»77). 

Next we find a recurrence for the sequence (3t ■ We start with the following 
lemma. 

Lemma 3. If U t n is preconfined, then U t+l rj(x) > 1 if and only if 

r)(x)eZ-(p t ,l3 t+1 }. 
Proof. By Lemma 2, since \J tq r\ is preconfined, we have 

r t < U t+1 7](x) < 1 + r t 
for all x G [0, 1], where rt = r{U t ri). Thus U t+1 rj(x) > 1 if and only if 

U t+1 rj(x) G Z+ [0,r t ). 
By (10), this condition is equivalent to 

ri(x) + pt +1 €Z + [0,r t ). 
Using the fact that fit+i = Pt + r t, this in turn is equivalent to 
Ti(x)€Z+[-p t+1 ,r t -p t+1 )=Z-(p t ,p t+1 ]. □ 

The essential information in a generalized chip configuration r/ is contained 
in the associated probability measure fx on [0, oo) given by 

for Borel sets A C [0, oo). If rj arises from a chip configuration a on K n , 
then [x([a, b]) is the proportion of vertices that have between an and bn 
chips. For y£l, write (, y /j, for the translated measure £, y u(A) = u(A — y). 
For yd, write n\y for the restricted measure n\y(A) = ll(A n Y). Then 
the update rule (8) takes the form 

up = e^Hci)) + e^-VliLooj). (ii) 

Now consider the measure v on R given by 

!/(>!) = ^ ^(m - A) = ^ £{x | m - n(x) G A}. (12) 

mgZ m£Z 

By Lemma 3, if U tr q is preconfined, then 

r(C/ t+1 r/) = u(Z - (A, A+i]) = KA>A+l]> 

hence 

A +2 = A+i + KA,A+i]- (13) 
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This gives a recurrence relating three consecutive terms of the sequence fit- 
Our next lemma simplifies this recurrence to one relating just two consecu- 
tive terms. 

Lemma 4. If r] is preconfined, then for all t > 

A+i = fifit), 

where f : [0, oo) — > [0, oo) is given by 

/(y) = A + 1/(0, y] (14) 

and v is given by (12). 

Proof. By Lemma 2, [7*7/ is preconfined for all t > 0, so from (13), 

t-i 



s=0 
t-l 

s=0 

= v(0,Pt}. 

Hence A+i = fifit)- n 
The function / appearing in Lemma 4 satisfies 
f(y + l) = f(y) + v(y,y + l] 
= f(y) + i 

for all y > 0. Thus it has a unique extension to a function / : M — > M 
satisfying f(y + 1) = f(y) + 1 for all y G M. 

Note that / is nondecr easing. If it is also continuous, then it has a well- 
defined Poincare rotation number [9, 13] 

t— >oo I 

which does not depend on y. Here /* denotes the t-fold iterate /*(y) = 

fif-Hy)), with / = id. 

Viewing the circle S 1 as K/Z, the map / descends to a circle map / : 
S 1 — > S 1 . The rotation number p(f) measures the asymptotic rate at which 
the sequence 

yJ(y)J(f(y)),--- 

winds around the circle. 

We define the activity of a generalized chip configuration r\ as the limit 

a(r\) = lim — . 

t— >oo t 

From Lemma 4 we see that if r\ is preconfined and / is continuous, then this 
limit exists and equals the rotation number of /. 

Lemma 5. If w is preconfined, then fit = f (0) for all t > 0. 
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Lemma 6. If w is preconfined and f is continuous, then a(rj) = /?(/). 

The conditions that rj is preconfined and / is continuous are most suc- 
cinctly expressed in terms of the associated probability measure: is sup- 
ported on [0, 2) and is non-atomic, that is, fJ-({y}) = for all y G [0, 2). 

To complete the commutative diagram (7), it remains to describe the 
dynamics U on lifts of circle maps. Call a map / : R — > M a monotone 
degree one lift if / is continuous, nondecreasing and satisfies 

f(y + 1) = f(y) + 1 (15) 

for all y G R. Let J- be the set of monotone degree-one lifts / : R — )■ R, 
equipped with the L°° topology. We define the update rule U : T — > T by 

Uf = R-f(o) °f°Rf(o) 

where for y G R, we write i? y for the translation x i— >• x + y. 

Let .M be the space of non-atomic probability measures [i on the inter- 
val [0, 2), equipped with the topology of weak-* convergence: fi n — > [i if and 
only if fj, n [0, y) /i[0, y) for all y G [0, 2). For /x G 7W, define <£(/i) : R -> R 
by 

= mC 1 ) 2 ) + ^ m ~ y ' m l 

where if 6 < a we define //(a, b] := —//[ft, a). 

Theorem 7. If [i £ M, then cp(n) G J 7 . TTie map ^ : Ai — > T is continuous 
and preserves the dynamics: U o <f> = <f> o U. 

Proof. (j>((j) is clearly nondecreasing in y, and since \i is non-atomic, is 
continuous. Moreover, for any y G R, 

0(A*)(y + !) - 0(^)(y) = X] -y- l ' m -y] = M R ) = 1 

so (j>(fj,) G J 7 . 

If /U n G and \i n — ?■ // G M, then the cumulative distribution functions 
Fn(y) = Mn[0> y) converge pointwise to F(y) := /x[0, y). Since F n and F are 
continuous and nondecreasing, this convergence is uniform in y. Since \i n 
and // are supported on the interval [0, 2), we have for y G [0, 1] 

0O«n)(y) = M 1 , 2 ) + - V, 1] + Mn(2 - y, 2] 
= 2F n (2)-F n (l-y)-F n (2-y). 

The right side converges uniformly in y to 2F(2) — F(l — y) — F(2 — y) = 
(f>(fj,)(y). Since (f>(fJL n ) and (f>(fJ,) are degree one lifts, we have 

sup \<j>{nn){y) - 4>(v)(y)\ = sup \(j)(nn)(y) - <j>(p)(y)\ -> o 

j/eR 2/6[0,l] 

as n — )• oo, which shows that <j) is continuous. 
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Let P = n[l,2) = 0(m)(O). Then 

{Uo4>)^){y) = {R_po<j>^)oR p ){y) 



^2 K m - y - 



Write n = no + fii, where no = A*l [0,1) an d Ml = A* I [1,2) - By (H)> we have 

c/ M =e%+e^"Vi 

hence 

Ufi[l, 2) = /xo[l - P, 2 - 0) + Ml [2 - /3, 3 - /3) 
= M [l-Al) + M2-/3,2) 

mez 

where in the last line we have used that fi is non-atomic and supported on 
[0,2). Moreover 



^2 u /u(m -y,m} = ^2 no(m-y- f3,m- (3} + 

mGZ m£Z 

+ ^^i(m-y+l-/3,m+l-/3] 

meZ 

= -y- P,™- p\. 

m£Z 

We conclude that 

(0 o = 2) + £ l//i(m - y, m] 

meZ 



so 6 o U = U o 6. □ 



One naturally wonders how to generalize the construction described in 
this section to chip-firing on graphs other than K n . A key step may involve 
identifying invariants of the dynamics. On K n , these invariants take a very 
simple form: by (10), for any two vertices v, w € [n], the difference 

U l a(v) - U t a(w) mod n 

does not depend on t. Analogous invariants for parallel chip-firing on the 
n x n torus are classified in [6], following the approach of [8]. 
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3. Devil's Staircase 

Let /, f n , g be monotone degree one lifts (15), and denote by /, f n ,g the 
corresponding circle maps S 1 — > S 1 . Write / < g if f(x) < g(x) for all 
x G M, and / < g if f(x) < g(x) for all x G K. We will make use of the 
following well-known properties of the rotation number. For their proofs, 
see, for example [9, 13]. 

• Monotonicity. If / < g, then p(f) < p(g). 

• Continuity. If sup \f n - f \ ->■ 0, then p(f n ) -)■ p(f). 

• Conjugation Invariance. If g is a homeomorphism, then p(go f o 

g- 1 ) = p[f)- 

• Instability of an irrational rotation number. If p{f) ^ Q, and 
fl < f < h, then p(/i) < p(f) < p(f 2 ). 

• Stability of a rational rotation number. If p{f) = p/q G Q, 
and f q Id : S 11 — > S 1 , then for sufficiently small e > 0, either 

p(g) = p/q whenever / < g < f + e, 

or 

= p/q whenever / — e < g < j '. 
Let o"2, o"3, . . . be a sequence of chip configurations, with a n defined on K n . 
Suppose a n is stable, i.e., 

0<a n (v)<n-l (16) 

for all v £ [n]. Moreover, suppose that there is a continuous function F : 
[0, 1] [0, 1], such that for all < x < 1 

-#{v G [n] I aJv) < nx} Fix) (17) 
n 

as n — t- oo. Since the left side is nondecreasing in x, this convergence is 
uniform in x. We define $ : M — >• M by 

= [ x ] -F(\x] -x). (18) 

Note that (16) and (17) force F(0) = and F(l) = 1; by compactness, F is 
uniformly continuous on [0, 1], and hence $ is uniformly continuous on K. 

The rescaled activity phase diagram of cr n is the function s n : [0, 1] — > [0, 1] 
defined by 

s n {y) = a(a n + ny). 
As n — > oo, the s n have a pointwise limit identified in our next result. 

Theorem 8. If (16) and (17) hold, then for each y G [0, 1] we have 

s n (y) -»■ s(y) :=p(RyO<s>) 

as n —> oo, where $ is given by (18), and R y {x) = x + y. 

In Proposition 10, below, we show that under an additional mild hy- 
pothesis, the limiting function s(y) is a devil's staircase. Examples of these 
staircases for different choices of F are shown in Figure 2. 
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Figure 2. The devil's staircase s(y), when (a) F(x) is given 
by (6); (b) F(x) = \fx for x G [0,1]; and (c) F(x) = x + 
2~ sin27rx. On the horizontal axis y runs from to 1, and on 
the vertical axis s(y) runs from to 1. 



To prepare for the proof of Theorem 8, let r) n ^ y be the generalized chip con- 
figuration associated to a n + ny, let fj, n>y be the associated measure on [0, 2), 
and let f n>y = 4>{n n ^ y ) be the corresponding circle map lift. 

Lemma 9. Let y £ [0,1]. If (16) and (17) hold, then as n — > oo 

SUP \f n ,y ~ $ ° Ry\ ~> 

uniformly in y. 

Proof. The measures /x H) o are nonatomic, and by (16) they are supported 
on [0, 1), so for x G [0, 1] we have 

~y) = Mn l2 /[1,2) + ^ - x + y, m] 

-meZ 

= Mn,o[l - y, 2 - y) + ^ fJti,o(m -x,m-y] 
= Mn,o[l - y, 1) + Mn,o(l ~x,l-y] 

= Hn,o[l ~ X, 1) 
= fnfi( x )- 

Since / n>2/ and f n $ o i? y are degree one lifts that agree on [— y, 1 — y], they 
agree everywhere. 

By (17) we have fj, n> Q — > fi, where n([a,b]) = F(b) — F(a). Since \i is 
supported on [0, 1), we have for x £ [0, 1] 

4>{n){x) = - x, 1] = 1 - F(l - x) = $(x). 
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Since $ and 4>(p) are degree one lifts that agree on [0, 1], they agree every- 
where. Hence 

SUP \fn iV - $ O Ry\ = SUP \(f nfi - $) O Ry\ 

= sup|/ n)0 - $| 

= sup |</>(/% )0 ) - 0(m)|- 

The right side does not depend on y, and tends to zero as n — > oo by the 
continuity of (f) (Theorem 7). □ 

Proof of Theorem 8. By Lemma 6, Lemma 9 and the continuity of the ro- 
tation number, 

Sn{y) = a{r] rh y) = p{f n ,y) ~> p($ O R y ) 

By the conjugation invariance of the rotation number, p(<&oR y ) = p(R y o§), 
which completes the proof. □ 

Write & y = i^of, and let $> y : S 1 — > S 1 be the corresponding circle 
map. We will call a function s : [0, 1] — > [0, 1] a devil's staircase if it is 
continuous, nondecreasing, nonconstant, and locally constant on an open 
dense set. Next we show that if 

(§ y ) q ^Id for all y € [0, 1) and all q G N, (19) 

then the limiting function s(y) in Theorem 8 is a devil's staircase. 

Proposition 10. The function s(y) = p{Q y ) continuous and nondecreasing 
in y. If z G [0,1] is irrational, then s^ 1 (z) is a point. Moreover, if (19) 
holds, then for every rational number p/q G [0, 1] the fiber s~ 1 (p/q) is an 
interval of positive length. 

Substantially similar results appear in [9, 13]; we include a proof here for 
the sake of completeness. 

Proof. The monotonicity of the rotation number implies that s is nonde- 
creasing. Since sup|$ y — $> y i\ = \y — y'\, the continuity of the rotation 
number implies s is continuous. 

Since 3>(m) = m for all m G Z, we have s(0) = and s(l) = 1. By 
the intermediate value theorem, s : [0, 1] — > [0, 1] is onto. If z = s(y) is 
irrational, and y\ < y < y%, then & yi < & y < & y2 , hence 

s{yi) < s(y) < s(y 2 ) 

by the instability of an irrational rotation number. It follows that s^ 1 (z) = 
{.</}• 

If s(y) = p/q G Q, then by the stability of a rational rotation number, 
since {& y ) q ^ Id, there exists an interval / of positive length (take either 
/ = [y — e, y] or I = [y, y + e] for small enough e) such that p(Q y >) = p/q for 
all y' G /. Hence I C s~ 1 {p/q). □ 
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Note that if y < y' , then f n ^ y < f n ,y', so s n (y) = p(f n ,y) is nondecreasing 
in y. By the continuity of s, it follows that the convergence in Theorem 8 is 
uniform in y. 

Our next result shows that in the interiors of the stairs, we in fact have 
s n (y) = s(y) for sufficiently large n. 

Proposition 11. Suppose that (16), (17) and (19) hold. Ifs^^p/q) = [a,b], 
then for any e > 

[a + e,b-e] C s~ l (p/q) 

for all sufficiently large n. 

Proof. Fix e > 0. By Lemma 9, for sufficiently large n we have 

sup \f nt y - $ o R y \ < e 

for all y G [0, 1]. Then f n;V > R- e o $ o R y , so for any y > a + e we have by 
the monotonicity and conjugation invariance of the rotation number 

Sn(y) = P(fn,y) > p(R-e ° $ ° Ry) 

= P{Ry-e o $) = s(y - e) > s(a) = p/q. 
Likewise f n ^ y < R e o $ o R y , hence for any y < b — e 

Sn(y) = P(fn,y) < p(Re o $ o R y ) 

= p(R y+e o $) = s (y + e) < a(6) = p/ g . □ 

4. Short Period Attractors 

In this section we explore the prevalance of parallel chip-firing states 
on K n with small period. 

Definition. A chip configuration a on K n is preconfined if it satisfies 

(i) o~(y) < 2n — 1 for all vertices v of K n . 
If, in addition, a satisfies 

(ii) max„ o~(v) — mm v o~(v) < n — 1 
then a is confined. 

Equivalently, a is preconfined (confined) if and only if the generalized chip 
configuration r/(a) is preconfined (confined) as defined in section 2. 

Recall (2) that for any chip configuration a on K n , the sequence (J7*cr)t>o 
is eventually periodic. Denote its transient length by to; that is 

t = min{t > | U t a = U 1 ' a for some t' > t}. 

We say that "v fires at time t" if ZPaiv) > n. 

Lemma 12. If a(o~) < 1, then U l a is confined for all t > t^. 
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Proof. If a(o~) < 1, then at each time step, some vertex does not fire. If a 
given vertex v fires at time t, then 

U t+1 a(y) < U l a{v). 

Since U t a(v) £ Z>o for all t, there is some time time T v at which v does 
not fire. By Lemma 2, we have U l a{v) < 2n - 1 for all t > T v , and U* a is 
confined for all t > max B T v . 

For any t > to we have iPa = U a for infinitely many values of t' , 
hence U l a is confined. □ 

For a chip configuration a on K n and a vertex v £ [n], let 

n t (fj, v) = #{0 < s < 1 1 C/V(w) > n} 

be the number of times v fires during the first t updates. During these 
updates, the vertex v emits a total of nut(a, v) chips and receives a total of 
at = Yl w u t{ (J i w ) chips, so that 

U^iv) - a(v) = a t - nu t {a,v). (20) 

An easy consequence that we will use repeatedly is the following. 

Lemma 13. A chip configuration a on K n satisfies U l a = a if and only if 

u t (a,v) = ut{a,w) (21) 

for all vertices v and w. 

Proof. If (21) holds, then ut(a, v) = ut/n for all v, so U t a(v) = a(v) by (20). 
Conversely, if £/V = a, then the left side of (20) vanishes, so ut(cr, v) = at/n 
for all vertices v. □ 

According to our next lemma, if a is confined, then ut(o~,v) and ut(o~,w) 
differ by at most one. 

Lemma 14. If a is confined, and o~(v) < o~(w), then for all t > 

u*(o-, v) < u t (a, w) < u t (a, v) + 1. 

Proof. Induct on t. If ut(a,w) = ut(a,v), then by (20) 

fjV('y) - C/V(io) = <r(v) - a(w) < 0. 

Thus if v fires at time t, then uu must also fire, hence 

u t+1 (a, v) < u t+ i(a,w) < u t+ i(a,v) + 1. (22) 

On the other hand, if Ut(cr,w) = Ut(cr,v) + 1, then since a is confined, we 
have by (20) 

U l a{v) - tfa^w) = n + a(v) - a{w) > 0. 
Thus if w fires at time t, then v must also fire, so once again (22) holds. □ 
Lemma 15. If a is confined, then U t o- = a if and only if n\at- 
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Proof. If U l a = a, then n\ott by Lemma 13. For the converse, write £ 
mm v ut(o~, v). By Lemma 14 we have 

u t {a,v) G {1,1+1} 

for all vertices v. If nlaj, then since 



D = l 



we must have ut(o~, v) = i for all v , so f/V = cr Lemma 13. □ 

Let a be a confined state on K n . By the pigeonhole principle, there exist 
times < s < t < n with 

a s = at (mod n). 

By Lemma 15 it follows that U s o~ = U l o~, so a has eventual period at most n. 

Our next result improves this bound a bit. Write m(a) for the eventual 
period of a, and 

u(a) = #{a(v)\v £ [n]} 
for the number of distinct heights in a. 

Proposition 16. For any chip configuration a on K n , 

m(o~) < v(o~). 

Proof. If a(o) = 1, then m(a) = 1. Since u{Ua) < u(a) and m{Ua) = 
m(o~), by Lemma 12 it suffices to prove the proposition for confined states a. 
Moreover, by relabeling the vertices we may assume that cr(l) > cr(2) > 
. . . > o~(n). Write the total activity at as 

a t = Qtn + R t 

for nonnegative integers Qt and Rt with Rt < n — 1. By Lemma 14, since 
ut(cr, v) = at, we must have 



u t (a,v) 



Qt + 1, v<R t 
Qt, v > R t , 



Moreover if R t > 0, then a{Rt) > cr(Rt + 1). Since the set 

V = {0} U {v G [n - l]\cr(v) > a(v + 1)} 

has cardinality v((j), by the pigeonhole principle there exist times < s < 
t < v{a) with R s = R t . Then 

ut(cr, v) - u s (a, v ) = Q t - Q s 

independent of v, hence U s a = U t a by Lemma 13. Thus m(a) < t — s < 
v{a). □ 
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Bitar [ ] conjectured that any parallel chip-firing configuration on a con- 
nected graph of n vertices has eventual period at most n. A counterexample 
was found by Kiwi et al. [11]. Anne Fey (personal communication) has found 
a counterexample on a regular graph. Proposition 16 shows that Bitar's con- 
jecture holds on the complete graph; it is also known to hold on trees [4] 
and on cycles [7]. It would be interesting to find other classes of graphs for 
which Bitar's conjecture holds. 

Next we relate the period of a chip configuration to its activity. We will 
use the fact that the rotation number of a circle map determines the periods 
of its periodic points: if / : R — > M. is a monotone degree one lift (15) with 
p(f) = p/q m lowest terms, then all periodic points of / : S* 1 — > S 1 have 
period q; see Proposition 4.3.8 and Exercise 4.3.5 of [9]. 

Note that if a chip configuration a on K n is periodic of period m, and 
a(ci) = p/q, then 

^ = V - (23) 
ran q 

and akm = ka m for all k £ N. 

Lemma 17. If a{a) = p/q and (p,q) = 1, then m(o~) = q. 

Proof. Since a(Ua) = a(o~) and m(Uo~) = m(o~), we may assume that a is 
periodic. Now if a(a) = 1, then m(a) = 1, so we may assume a(a) < 1. In 
particular, a is confined by Lemma 12. Write m = m(a). By Lemma 13 
and (23) 

a m pm 

Um(a,v) = = 

n q 

for all vertices v. Hence q\pm and hence q\m. 

Let / be the lift of the circle map / associated to a. By Lemma 6 we 
have p(f) = a(a), so all periodic points of / have period q. By Lemma 5 we 
have 

fqm{0) = ^rn = qOrn =pm£ ^ 

n n 

so is a periodic point of /. Thus f q (0) £ Z. Now by Lemma 5 

a q = n/«(0) 

hence U q a = a by Lemma 15. Thus m\q. □ 

Given 1 < p < q < n with (p, q) = 1 and p/q < 1/2, one can check that 
the chip configuration a on K n given by 



a(v) 



v+p—1, v<q — l—p 

v + n + p — q—1, q—p<v<q 
n + p — 1, v > q. 



has activity a(o~) = p/q. For a similar construction on more general graphs 
in the case p = 1, see [14, Prop. 6.8]. By Lemma 17, m(o~) = q. So for every 
integer q = 1, . . . , n there exists a chip configuration on K n of period q. 
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Despite the existence of states of period as large as n, states of smaller 
period are in some sense more prevalent. One way to capture this is the 
following. 

Theorem 18. If 02, 03, . . . is a sequence of chip configurations satisfying 
(16), (17) and (19), then for each q E N there is a constant c = c q > such 
that for all sufficiently large n, at least cn of the states {a n + k}^ =0 have 
eventual period q. 

Proof. Fix e > 0, and choose p < q with (p,q) = 1. By Proposition 10 we 
have s~ 1 (p/q) = [a p , b p ] for some a p < b p . By Proposition 11, for sufficiently 
large n we have 

a(a n + [ny\) = s n {y) =p/q 
for all y S [a p + e/2q, b p — e/2q]. Thus all states a n + k with 

a p+ ^jn<k<(b p --^n 

have activity p/q and hence, by Lemma 17, eventual period q. We can 
therefore take 

p:(p,q) = l 

for any e > 0. □ 

The rest of this section is devoted to proving the existence of a period 2 
window: any chip configuration on K n with total number of chips strictly 
between n 2 — n and n 2 has eventual period 2. 

The following lemma is a special case of [ I , Prop. 6.2]; we include a proof 
for completeness. Write r(cr) for the number of unstable vertices in a. 

Lemma 19. // a and r are chip configurations on K n with o~{v) + t{v) = 
2n — 1 for all v, then a(a) + a(r) = 1. 

Proof. If o~(v) + t(v) = 2n — 1 for all v, then for each vertex v, exactly one 
of a{v ), t{v) is > n. Hence r{a) + r(r) = n, and 

Ua{y) + Ut(v) = a(v) + r(cr) + t(v) + r(r) — n 
= 2n-\ 

for all vertices v. Inducting on t, we obtain 

r(E/V) + r(C7*r) = n 

for all t > 0, and hence a(o~) + a(r) = 1 by (3). □ 

Given a chip configuration a on K n , for j = 1, . . . , n we define conjugate 
configurations 

/ \ J 0-0) + j - n, v < j 
[a{v)+j, v>j. 

Note that c n a = a. 
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Lemma 20. Let a be a chip configuration on K n , and fix j G [n]. For 
all t > 0, we have for v < j 

u t (a,v) - 1 < Ut(c?o-,v) < u t (o-,v), 

while for v > j 

uticr, v) < u t (c?o-, v) < ut(a, v) + 1. 
Proof. Induct on t. By (20) 

[/Vu(w) = cPa(v) + a t> j - nu t (c? a,v) (24) 

where 

n 

a t ,j = Ut{c>o-,w). 

w=l 

By the inductive hypothesis, 

- j < a t j - a t ,o <n- j. (25) 

Fix a vertex v, and write 

b = Ut(c i o-,v) - u t (a,v). 

By the inductive hypothesis we have b G {—1, 0, 1}. From (24) we have 

t/V 'a{v) - L/V(v) = a t j - a t ,o + cP a{v) - a(v) - nb (26) 

Case 1: v < j and b = — 1. By (25), the right side of (26) is > 0, so if v 
is unstable in U t a, then it is also unstable in IPcPa. Thus 

u t +i(o-,v) - 1 < ii i+ i(c J o-,'t;) < u t+1 (a,v), (27) 

completing the inductive step. 

Case 2: v < j and 6 = 0. By (25), the right side of (26) is < 0, so if v 
is unstable in U t c :, a, then it is also unstable in U t a. Thus once again (27) 
holds. 

Case 3: v > j and 6 = 0. By (25), the right side of (26) is > 0, so 

u t+1 (a,v) < u t+ i{c>a,v) < u t +i(a,v) + 1, (28) 

completing the inductive step. 

Case 4: v > j and 6 = 1. By (25), the right side of (26) is < 0, so once 
again (28) holds. □ 

Corollary 21. For any chip configuration a on K n and any j G [n], 

a(c'a) = a(a). 

It turns out that the circle maps corresponding to a and cPa are conjugate 
to one another by a rotation. This gives an alternative proof of the corollary, 
in the case when both a and c 3 o~ are preconfined. 

Lemma 22. Let a be a chip configuration on K n . Lfu2(cr,v) > 1 for all v, 
then v,2t(cr, v ) > t for all v and all t > 1. 
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Proof. Induct on t. By the inductive hypothesis, a^t > nt. Fix a vertex v, 
and suppose that U2t(o~,v) = t. Since U2(cr,v) > 1, either a(v) > n or 
Ua(v) > n. In the former case, by (20) 

U 2t a(v) - a(v) = a 2t - nt > 0, 

so f fires at time 2t. Hence v,2t+i(o~, v) > t + 1. Summing over v yields 

&2t+i > £*i(t + 1) + (n — «i)t = nt + ai. 

In the latter case, if v does not fire at time 2t, we have by (20) 

U 2t+1 a(v ) - Ua(v) = a 2t +i - «i - nt > 0, 

so f fires at time 2t + 1. Hence U2t+2(o~, v) > t + 1 for all v. □ 

Given a chip configuration a on i^ n , write 

£(a) = min{cr(l), . . . , cr(n)} 

and 

r(a) = #{v G [n] : cr(i?) > n}. 

Write 

n 
v=l 

for the total number of chips in the system. 

Theorem 23. Every chip configuration a on K n with n 2 — n < \o~\ < n 2 
has eventual period 2. 

Proof. If \o~\ < n 2 , then a(o~) < 1, so by Lemma 12 we may assume a is 
confined. Moreover, by relabeling the vertices we may assume that o~(l) > 
<t(2) > . . . > a(n). In particular, for any vertex v, since o~(w) > o~(v) — n + 1 
for w = v + 1, . . . , n, we have 

\a\ > va(v) + (n — v)(a(v) — n + 1) 
= na(v) — (n — v)(n — 1). 

Since \o~\ < n 2 , we obtain 

, , n 2 + (n — v )(n — 1) 

a(v) < — ^ <2n-v. 29) 

n 

For a fixed vertex v, as j ranges over [n], the quantity c^a(v) takes on 
each of the values 

a(v) + v — n, . . . , cr('w) + t> — 1 
exactly once. By (29), at least o~(v) + v — n of these values are > n, hence 

r(Jo-) > "£(a(v) + v-n) = \a\- l^L^. ( 30 ) 

j=l v=l 
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Since £(c 3 a) = a(j) + j — n, we have 

± e{ ^ ) = M .^-21. (31) 

Summing equations (30) and (31), and using |a| > n 2 — n, we obtain 

n 

^(^(cV) + r(cV)) > n 2 — n. 

3=1 

Since each term in the sum on the left is a nonnegative integer, we must 
have 

£{c?o) + r(cPcr) > n. 

for some j G [n] . Thus every vertex v fires at least once during the first two 
updates of cPa. From Corollary 21 and Lemma 22 we obtain 

a(a) = a(cPa) > -. 

Finally, the chip configuration t{v) := 2n — 1 — <j(v) also satisfies n 2 — 
n < \t\ < n 2 , so a(r) > g- By Lemma 19 we have a(a) + a(r) = 1, so 
a(a) = a(r) = \. From Lemma 17 we conclude that m(a) = 2. □ 
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